Resonance behavior of the defect-induced Raman mode of single-chirahty enriched 

carbon nanotubes 
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We present a resonance Raman study of the disorder-induced D mode in a sample highly enriched 
with semiconducting (9,7) single-walled carbon nanotubes in the excitation energy range of 1.49 — 
2.05 eV. The intensity of the D mode shows a resonance behavior near the optical transition of the 
(9,7) tube. The well-known dispersion of the _D-mode frequency, on the other hand, is not observed 
at the resonance, but only above a certain excitation energy. We explain our results by numerical 
simulations of the _D-mode spectra. 

PACS numbers: 



I. INTRODUCTION 



Raman spectroscopy is a powerful characterization tool 
for carbon nanotubes and graphene, due to its high sen- 
sitivity to sp^ -hybridized carbon [l|-|3|. In particular, the 
disorder-induced D mode and the related two-phonon 2D 
{G') mode are indicative of sp^-hybridized carbon, as 
they stem from a breathing-like vibration of the carbon 
hexagons Q in a double-resonant Raman process (dRR) 
Q . Because of their origin the D and 2D modes are most 
sensitive to the electronic band structure as well. This, 
among others, leads to the spectroscopic distinction be- 
tween single- layer and few-layer graphene Q. 
The disorder-induced D mode, in addition, depends on 
the concentration of structural defects in carbon nan- 
otubes and other graphitic materials 0|; its intensity 
is therefore taken as evidence for, e.g., successful cova- 
lent functionalization d, @ ■ In single- walled carbon nan- 
otubes (SWCNTs), however, the intensity and resonance 
behavior of the D mode have not yet been fully under- 
stood [Tol . [ll| . This is partly due to the fact that, initially, 
chiral-index defined SWCNT samples were not available. 
Moreover, most experiments on the D mode of carbon 
nanotubes investigated the dispersion of the f-mode fre- 
quency with excitation energy. They were done either in 
nanotube ensembles, using a broad range of excitation 
energies [l^ - [T5| . or on individual SWCNTs, using only 
a small number of excitation energies and assuming that 
only those SWCNTs contribute to the D mode for which 
the excitation energy is in resonance with an optical tran- 
sition 

EMI. 

Numerical simulations, on the other hand, 
focused on the D-mode dispersion of both, individual 
nanotubes [l3| and nanotube ensembles [l5| . showing an 
intensity decrease of the D mode away from the optical 
transition of the nanotube. Ref. [l5| emphasized the in- 
tensity enhancement of the D mode directly at optical 
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resonances (enhanced-dRR scattering), attributing the 
D-Tnode dispersion in nanotube ensembles to excitation- 
energy dependent contributions of different nanotubes 
with different D-mode frequencies. However, excitation- 
energy dependent Raman experiments following the D- 
mode dispersion and intensity on chiral-index identified 
carbon nanotubes have not been reported to the best 
of our knowledge. On the other hand, knowledge about 
nanotube phonon dispersion [l9l [20|. optical transitions 
(2l| - [23t . and chiral-index sorting and assignment (234271 
has improved significantly over recent years. 
Here we present a resonance Raman study of the D mode 
in a sample highly enriched with (9,7) SWCNTs in the 
excitation energy range of 1.49 — 2.05 eV. The intensity 
of the D mode in our sample depends strongly on excita- 
tion energy and shows a resonance behavior near the op- 
tical transition of the (9,7) SWCNT. On the other hand, 
the well-known dispersion of the i?-mode frequency with 
excitation energy starts only above a certain excitation 
energy, i.e. away from the resonance. Numerical simula- 
tions of the D-mode spectra of the (9,7) SWCNT agree 
with both observations, supporting the interpretation of 
an enhanced, only weakly dispersive double-resonant Ra- 
man scattering at excitation energies close to the optical 
transitions in SWCNTs. 



II. EXPERIMENTAL SETUP 

SWCNTs produced by pulsed-laser vaporization were 
wrapped with the polymer Poly(9,9-di-n-octylfuorenyl- 
2,7-diyl) (POF) in Toluol. After ultracentrifugation 
mainly individual (9,7) SWCNTs were dissolved in this 
solution [2^. After a time of rebundling, the sample was 
centrifuged again to separate the (9,7) SWCNTs from 
excessive POF. This was necessary to reduce the strong 
Raman signal of the POF in the frequency region of the 
D mode. The obtained solution was drop casted onto a 
Si/SiOa (200 nm SiOa) substrate. 

Raman measurements were performed with a microscope 
setup in backscattering geometry. A titanium-sapphire 
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Figure 1: (Color online) RBM resonance profiles of sample 
enriched with (9,7) SWCNTs with chiral-index assignment. 
The inset shows the Raman spectra of the G~ modes of the 
(9,7) and (9,8) SWCNTs at 1.55 eV excitation energy. 

(1.49 - 1.77 eV) and a dye laser (1.79 - 2.05eV) were 
used for excitation with 1 mW laser power. The spectra 
were collected using a Dilor XY800 triple monochroma- 
tor. They were calibrated with a neon lamp; the intensity 
was normalized to CaF2 for the radial-breathing mode 
and to diamond for the D mode. All measurements were 
performed at the same spot of the sample. 



III. EXPERIMENTS 

A. Chiral-index assignment 

Figure [1] shows the resonance profiles of the radial- 
breathing mode (RBM) for the chiral-index assignment 
[28| . The profiles with the highest intensities can be as- 
signed to the second optical transition {E22) of semicon- 
ducting (9,7) and (9,8) SWCNTs. We found that the 
RBM of the (9,7) SWCNTs consists of two contributions 
with slightly different frequencies (Aw = 2.5 cm^^) and 
resonance energies {E' = 1.5405 eV, E" = 1.5336 eV, 
AE = 6.9 ± 0.8 meV). Assigning these resonance profiles 
to two different chiral indices seems not suitable, because 
of the isolated position of the (9,7) SWCNT with respect 
to other SWCNTs in the Kataura plot [Hi. A splitting of 
the RBM resonance profile due to incoming and outgo- 
ing resonance would correspond to an energy difference 
of approximately 27meV, which is much higher than in 
our measurements. Therefore, the splitting cannot be 
explained up to know in detail but might be due to the 
enrichment process. We confirm our chiral-index assign- 
ment with help of the G~ mode (see Table |l| [Hi. The 
Raman spectrum of the G~ modes at 1.55 eV excitation 
energy is shown in the inset of Fig. [T] 
Some weaker RBMs of other SWCNTs (Fig.[Tl open sym- 



SWCNT G--mode frequency (cm-^) 

this work exp. Ref. [26] theory Ref. [26| 

(9.7) 1556.7±0.8 1557.4 ± 1 1558.8 

(9.8) 1564.4±0.7 - 1561.3 



Table L G"-mode frequencies of the (9,7) and (9,8) SWCNTs 
compared with Ref. [2a |. 



bols) are observed above the excitation energy of 1.6 eV, 
however, their intensities were at least by a factor of 27 
lower than for the (9,7) SWCNT and by a factor of 3 
lower than for the (9,8) SWCNT. 



B. Dispersion and resonance behavior of the 
defect-induced Raman mode 

Figure [5] shows the Raman spectra of the D mode for 
excitation energies between 1.49 and 2.05 eV, where the 
maximum amplitude was normalized to one. In these 
spectra, the D mode consists of at least two main contri- 
butions. They are indicated by two colored stripes and 
denoted with the signs + and x . For high excitation en- 
ergies, the Raman signal of POF appears in the spectra 
due to a decrease of the D-mode intensity. The spectra 
were fitted with two Lorentzians, after the background 
signal of the POF polymer was subtracted. 
Figure [3] shows the measured I?-mode intensity (peak 
area) of both contributions (+, x) as function of exci- 
tation energy on logarithmic scale. Both contributions 
show a decrease in intensity by two orders of magnitude 
for the excitation energy ranging from 1.49 to 2.05 eV. 
The dispersion of the D-vaoAe frequency with excitation 
energy is shown in Fig. [4l The frequencies are repre- 
sented by + and x signs, the intensities of the D mode 
contributions are indicated by the radii of the circles. 
The strongest D-raode intensity is observed, when the 
RBM of the (9,7) and (9,8) SWCNTs is in resonance. 
Unexpectedly, the well-known dispersion behavior of the 
D mode [l^ is only observable above a certain excitation 
energy (1.7 cV for the SWCNTs in our sample). Here, 
the D-xnode intensity is strongly decreased, see Fig. [3l 
the dispersions of the two contributions can be linearly 
fitted to 81cm-7eV and 43cm-VeV (Fig.g]). 
We explain our results tentatively by attributing the 
weakly dispersive, high-intensity appearance of the 
D mode (Fig. [H excitation energy between 1.49 and 
1.7 cV) to an enhanced-dRR scattering very close to the 
optical transition. For the dispersive, low-intensity part 
of the D mode (Fig. |4l excitation energy above 1.7eV), 
we are above the resonance window of the enhanced-dRR 
scattering, and therefore the intensity decreases. At the 
same time, the electron bands of SWCNTs exhibit the 
nearly linear dispersion derived from the nearly linear 
bands close to the if -points of graphite and graphene. 
To verify this interpretation, we performed calculation of 
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Figure 2: (Color online) Raman spectra of the D mode of the 
sample enriched with (9,7) SWCNTs for excitation energies 
between 1.49 and 2.05 eV (top to bottom). Amplitudes of the 
spectra were normalized to 1; the spectra are plotted verti- 
cally offset. Colored stripes and signs (+, x) indicate two 
main contributions of the D mode. For high excitation ener- 
gies, the Raman signal of POF (bottom spectrum) overlaps 
with the weak i3-mode intensity. 
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Figure 3: (Color online) Resonance behavior of the D mode 
of the (9,7) enriched sample for excitation between 1.49 and 
2.05 eV. The signs + and x denote the two contributions 
marked in Fig. [3 



1330 - 



^ 1320 

'e 

^ 1310 
c 

CO 

E 1300 

CO 

q: 

1290 



_^ 1 1 1 1 1 1 1 1 1 1 1 1 1 


1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 








X — 

_!^1 cm'VeV - 




^X"" + 




+ 42 cm"VeV- 




: 


1^22 


1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 


1.5 1.6 1.7 


1.8 1.9 2.0 



1280 



Excitation energy (eV) 

Figure 4: (Color online) Dispersion of the two contributions 
of the D mode, denoted by the signs + and x (Fig. [2]). The 
radii of the circles indicate the D-mode intensities. The ar- 
row indicates the E22 transition energy of the (9,7) SWCNT 
determined from the RBM resonance profile (Fig. [T|). 



the resonance behavior of the D mode for the (9,7) and 
the (9,8) SWCNT. 



IV. SIMULATION 



Simulation details 



The D mode is calculated by a numerical simulation of 
the dRR scattering. The differential cross section is pro- 
portional to |/Cp [30|. The matrix element /C for the scat- 



tering with one phonon and one defect is given by |l3| : 

a,b,c 

Me~rMe-defMe-phMe-r 

[E-Ef,- n)iE - El - hujj,^ - i-i){E - - h^j,t, - n) 

Me-rM^-phMe-defMe^r 



{E-El- i-i){E - El - i-i)(E - E^ - hwph - i-y) 

The first term describes the scattering process by first 
scattering the electron by a phonon and then the hole 
by a defect; the second term describes the scattering 
process by first scattering the electron by a defect and 
then the hole by a phonon. E is the excitation energy. 
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E'Kx = a,b,c) are the energy differences between ex- 
cited and ground state of the electronic bands, hujph is 
the phonon energy and 7 is a empirical broadening fac- 
tor taking into account the lifetime of the intermediate 
electronic states. The different A^'s are the matrix ele- 
ments for the interaction between electron and radiation 
(Me-r), phonon (Me-ph), or defect (Me-def)- Each M 
was set constant for the simulation; 7 was set to 40 meV. 
The electronic bands were calculated from a third nearest 
neighbor tight-binding approximation (3l] |. The phonon 
dispersion relations were calculated with a symmetry- 
based force constants approach adapted to the exper- 
imental phonon dispersion of graphite [2^ [32| with the 
program POLSym In chiral SWCNTs, both the lon- 
gitudinal optical (LO) and the transverse optical (TO) 
phonons are Raman active at the T point, in contrast 
to zigzag or armchair SWCNTs Inside the Brillouin 
zone, they lose their defined longitudinal or transverse 
character; therefore they will be called here LO- and TO- 
derived phonon branch. 

Figure [5] shows the calculated electron and phonon dis- 
persion of the (9,7) SWCNT in units of helical quantum 
numbers [l|, HJ- The transition energies En and £^22 
in the electron dispersion, and the LO- and TO-derived 
phonon branches in the phonon dispersion are denoted. 
In our simulation. En = 0.67cV and £'22 = 1.28eV [sst . 



B. Calculated dispersion and resonance behavior of 
the D mode 

We performed simulations with excitation energies in 
the range of both optical transitions En and i?22- There- 
fore, up to four scattering processes for the electrons are 
in principle possible: En o En, E22 ^ E22, En ^ E22 
and £^22 ^ En. A sketch of these processes for the 
lowest possible excitation energy is shown in Fig. [6l[a). 
Including the two phonon branches and the four scatter- 
ing processes, eight different processes contribute to our 
D-mode simulation [Fig. [6l^b)| . 

In contrast to Ref. [13j, where we assumed that the dis- 
order only breaks the pure translational symmetry and 
not, to first approximation, the rotational symmetry of 
the nanotube, we now allow the entire phonon branch 
with helical quantum number to = 0. This corresponds 
in linear quantum numbers to a change in band index to, 
which would not be allowed if the rotational symmetry 
was conserved. 

To analyze differences between the various contributions 
of the simulated Raman spectra, we first calculated each 
scattering process separately, thereby neglecting possi- 
ble interference effects. Each spectrum of the differ- 
ent contributions was phenomenologically fitted with two 
Lorentzians (inset Fig. [7]) considering the shorter and 
longer phonon wave vector from the scattering processes 
at a given energy. 

Fig. [7] shows the fitted data of the simulated D mode for 
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Figure 5: (Color online) Calculated electron and phonon dis- 
persion of the (9,7) SWCNT. The optical transitions £11 and 
£22, and the LO- and TO derived phonon branches are indi- 
cated. The wave vector is given in units of helical quantum 
numbers with q = 386, n = 1 and lattice translation 

vector a = 59.22 A. 



each of these eight different contributions. Different col- 
ors denote the different scattering processes (color code 
as in Fig. [6]) and the shapes of the symbols denote the 
used phonon branches. The symbol size indicates the in- 
tensity obtained from the fit. 

From these results, some general observations can be 
made. First, all scattering processes show an enhanced 
intensity near the corresponding optical transitions (£'11, 
£^22)- For the symmetric scattering within the same opti- 
cal transition (£'11 ^ En and £^22 ^ £'22), the D-mode 
intensity of the outgoing resonance (w 160 meV above the 
optical transition) is always stronger than for the incom- 
ing resonance. The reason is that for the incoming reso- 
nance, the electron scattered with a phonon will always 
be below the optical transition (i.e. band extrema with 
high density of states in our simulation) , whereas for the 
outgoing resonance, scattering with a phonon matches 
the optical transition. Second, the D mode becomes dis- 
persive at energies above the outgoing resonance, as ob- 
served in experiments. 

A peculiarity in the dispersion is the scattering process 
for the shorter wave vector of the scattering £^22 E22 
involving the TO-derived branch (Fig. [71 gray circles, 
lower branch). Here, an unusual negative dispersion can 
be observed. This is due to the fact that the phonon wave 
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Figure 6: (Color online) (a) Electron dispersion from Fig. [S] 
zoomed to the energy region of the optical transitions £11 and 
E22- The upper panel shows a sketch for the '"symmetric"' 
scattering {En -f^ En, E22 E22) of electron and hole; the 
lower panel for the antisymmetric scattering ( En •<->■ £22, 
i?22 •<->■ Ell), (b) The corresponding phonons for the scatter- 
ing processes in (a) are drawn into the LO-derived (upper) 
and the TO-derived branch (lower panel). 



vector q of the lowest possible optical transition does not 
correspond to the lowest frequency of the TO-derived 
phonon branch [Fig. [S]Jb), lower panel, gray arrow]. For 
smaller q vectors (corresponding to higher excitation en- 
ergy) , therefore the frequency of the D mode first drops 
before it rises again after passing the minimum of the 
phonon dispersion. 

A full simulation of the D-mode intensity of the (9,7) 
SWCNT with an overlap of the eight different scatter- 
ing processes, necessarily including all interference effects 
(3g, is shown in Fig. [5]Ja). An excitation-energy range 
similar to the experiment was chosen. 
Similar to the experiment, the simulated D mode shows 
a strong intensity in a small energy range, followed by an 
intensity decrease and a starting dispersive behavior for 
higher energies. A comparison with the calculated data 
from Fig. [7] shows that the highest intensity belongs to 
the scattering E22 ^ E22 with the LO-derived phonon 
branch (gray squares). 

To deduce the dispersion of the full simulation, the Ra- 
man spectra were phenonienologically fitted with five to 
six Lorentzians [inset of Fig. [5{b)]. The result is shown 
in Fig. [H^b). Here the radii of the circles correspond to 
the D-mode intensities. 

The inset of Fig. [SJ^b) shows that in the frequency range 
of the experimental data (1280 — 1330cm~^), the simu- 
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Figure 7: (Color online) Calculation of the different contribu- 
tions of D mode for the (9,7) SWCNT. Different colors denote 
the scattering path; the shapes of the symbols denote the used 
phonon branch. The symbol size indicates the intensity. The 
inset shows exemplary a fit of a simulated spectrum with two 
Lorentzians. 



lated Raman spectra consist of two main contributions. 
These two main contributions are due to an overlap of 
the different scattering processes. For these two main 
contributions a dispersion of 38 and 56cm~^/cV (solid 
lines) is obtained. 

For our calculation of the D-mode resonance of the (9,8) 
SWCNT (not shown) we found a similar behavior as de- 
scribed above. For that reason and due to the much 
weaker RBM intensity of the (9,8) SWCNT in compar- 
isons to the (9,7) SWCNT we neglected this contribution 
of the simulated D mode in the following comparison be- 
tween calculation and experimental data. 



V. COMPARISON OF CALCULATIONS WITH 
EXPERIMENTAL DATA 

In Fig. [21 the experimental data (Fig. ^ are plotted 
into the intensity color plot of our simulations [Fig.[8lja)]. 
In our simulations, the energies of the optical transitions 
were slightly underestimated. For better comparison, the 
simulations were shifted by -1-0.28 eV to bring the theo- 
retical optical transition energies to the same energy as 
in experiment. 

In general, the agreement between calculation and ex- 
periment is quite good as discussed above: In both cases 
an enhanced D-mode intensity is visible, as well as the 
non-dispersive parts have similar frequencies. Further, 
at higher excitation energies, we observe a similar dis- 
persion of the D mode. 

The two most obvious differences between simulation and 
experiment are: First, the D-mode intensity in the sim- 
ulation decreases only linearly by a factor of two and not 
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Figure 8: (Color online) (a) Simulated D-mode intensity as 
a color plot with consideration of all interference effects, (b) 
Dispersion of tlie simulated Raman spectra, deduced by fit- 
ting with Lorentzians (inset). The radii of the circles corre- 
spond to the simulated D-mode intensity. 



exponentially by two orders of magnitude; second, the 
highest D-mode intensity in the simulations is shifted 
to higher excitation energy. This shift is caused by the 
strong outgoing resonance in the simulation, which we do 
not observe in the experiment. Both differences might be 
due to neglecting the wave-vector dependence of matrix 
elements and excitonic effects in our simulations. Also 
an intensity decrease of the outgoing resonance due to 
the non-Condon effect (recently observed for the G mode 
(37| ) might play a role for the D-mode intensity and was 
not taken into account in our simulation. 

A comparison of the Raman spectra between experi- 
ment and simulation is shown in Fig. 1101 For excitation 
energies near the optical transition £22, a small differ- 
ence between experiment and simulation is seen. This 
frequency difference might be caused by inaccuracies in 
the calculated phonon dispersion, in particular near the 
wave vectors derivated from the K point in graphene, 
which contribute most at excitation directly at the op- 
tical transition. Since also the different contributions of 
the D mode have slightly different frequencies (Fig. [7]), 
the observed frequency deviations might also be caused 
by an overestimation of one of these contributions, again 
because of neglecting the matrix elements. 
Furthermore, an additional peak appears for higher ex- 
citation energies in the range of low frequencies in our 





1350 




1340 




1330 




lozO 


'c 


1 1 U 


0^ 






1300 




1290 








1280 


TO 




E 


1270 








1260 




1250 




1240 




1230 




1220 



experimental 




■ 10 

8 

6 

4 
,2 
'0 



2.4 



1.5 1.6 1.7 1.8 1.9 2.0 2.1 2.2 2.3 
Excitation energy (eV) 



Figure 9: (Color online) Comparison of dispersion behavior 
between simulations (color intensity plot) and experimental 
data (symbols). 
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Figure 10: (Color online) Comparison of calculated (dashed 
lines) and experimental (solid lines) Raman spectra. Ampli- 
tudes of the spectra were normalized to 1; the spectra are 
plotted vertically offset. 



simulations (Fig. [101 at w 1250 cm^^). This contribution 
corresponds to a scattering processes involving the TO- 
derived phonon branch (Fig. [71 circles) . Due to the anti- 
crossing of the phonon bands near the K point, the fre- 
quencies of the TO-derived phonon branch corresponds 
to the LO-phonon branch at the F-point for the (9,7) 
SWCNT. For F-point scattering (RBM and the G mode) 
it was shown in calculations that for SWCNTs with a 
large chiral angle [as for the (9,7) nanotube], the Ra- 
man intensity related to the LO-phonon branch is much 
weaker than for the TO-phonon branch There- 
fore we think that the Raman intensity for scattering 
involving the TO-derived phonon branch at the K point 
(Fig. [TUl at w 1250 cm"-'^) might be much smaller than 
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in our simulations; as we observe in our experiments. 
The unexpected high dispersion of 81cm^^/cV in the 
experimental data (Fig. 3]) might be explained with a 
jump to contributions from another scattering process 
as seen in the calculations in Fig. [St^b), indicated with 
a dashed line. Another explanation might be that due 
to the exponential decrease of the D mode the still re- 
maining metallic nanotubes contribute to the D mode for 
excitation energies above 1.9 eV. This would lead to an 
additional contribution to the D mode at higher frequen- 
cies [llj and therefore to an overestimated dispersion. 

VI. CONCLUSION 

In conclusion, we showed the resonance behavior of the 
Raman D mode in a sample enriched with (9,7) SWC- 
NTs, thereby following the D mode of very few, chiral- 
index identified nanotubes across an optical transition. 
The D mode intensity is strongest at the optical transi- 
tion of the (9,7) SWCNT; the D-mode frequency, how- 
ever, shows almost no dispersion with excitation energy 
near the resonance. The well-known dispersion of the D- 



mode frequency is observed only above resonance, where 
the intensity has dropped by approximately one order of 
magnitude. These results are understood based on nu- 
merical simulations of the £)-mode spectra, although the 
intensity drop observed in experiment is even stronger 
than in the simulations. Further calculations therefore 
should include excitonic effects. 

Based on these observations, we expect that the D-mode 
dispersion measured in nanotube ensembles not just re- 
flects the electron and phonon bands of a given nanotube 
but also a diameter and chiral-angle dependence of the 
Z)-mode frequency [Tsl 

Our results further show that the I?-mode intensity does 
not exclusively depend on the defect density of the nan- 
otubes but also on the excitation energy. This should 
be considered together with the excitation-energy depen- 
dence of the G-mode intensity [s^ when taking the D / G 
intensity ratio as a measure for the defect density. 
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